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Abstract. 

We report on the quantification of entanglement by means of entanglement measures on a 
four- and a six- qubit cluster state realized by using photons entangled both in polarization 
and linear momentum. This paper also addresses the question of the scaling of entanglement 
bounds from incomplete tomographic information on the density matrix under realistic 
experimental conditions. 



PACS numbers: 03.65.Ud, 03.67.Bg, 42.50.Ex 
Keywords: graph state experiment, entanglement measures 



Optimal verification of entanglement in a photonic cluster state experiment 



2 



1. Introduction 

Experiments in Quantum Information Science (QIS) rely heavily on multipartite entangled 
quantum states. Cluster states [Ij, or more generally graph states, are a particular class of 
multipartite states that offer a diversity of applications in QIS, ranging from measurement- 
based quantum computation and error-correction codes to nonlocality tests. Due to the 
importance of graph states, a considerable experimental effort has been made to realize them 
using photons [21 El IH El El ^ and cold atoms [8j . Proposals for trapped ions are also pursued 

In this paper we characterize the four- and six-qubit cluster states realized in [6l [71 [12] 
in terms of fidelity, purity and entanglement. In particular, we quantify the amount of 
experimentally generated entanglement using entanglement measures [13j. Cluster states are 
uniquely defined by a set of correlation operators which generate a group called the stabilizer. 
An interesting question this paper addresses is how bounds based on measurement results of 
the generator of the stabilizer alone scale with increasing system size under realistic conditions 
using the results developed in Refs. |T5l [TBj. We compare the optimal bounds based on 
such measurements for the fidelity, purity and the robustness of entanglement [TTl [181 HH] as 
well as the relative entropy of entanglement [20j with the density matrix obtained from all 
stabilizers in order to answer this question. 

2. Experimental Set-Up 

Cluster states are particular multiqubit entangled states associated to a graph. In the following 
we denote the Pauli spin matrices acting on the Hilbert space of qubit g by X^, Y^, and Zq. 
Given a lattice with n vertices and L links, a n-qubit cluster state can be defined by associating 
a qubit in the superposition state |+) = :^(|0) + |1)) to each vertex and a control-Z gate 
CZab = |0)a(0| (8) Ife + |l)a(l| ® ^6 to cach link between vertices a and b. In an equivalent 
way, the cluster state is defined as the unique eigenvector with positive eigenvalues of the 
n generators ga defined as ga = Xa YibeJVa where A/'a is the set of neighbouring vertices 
linked with a. The set of operators {ga} generate an Abelian group called the stabilizer S 
of the underlying graph. Note that eigenstates of the generators with negative eigenvalues 




Figure 1. left) Source of polarization-path hyper entangled state, right) Measurement setup 
for the path DOF. 
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9a\{i}) = ("l)^""!!^}) with \{i}) = 1^1, (the graph state basis states) are also referred to 

as graph states in the hterature. These states are equivalent to each other up to single qubit 
unitary transformations and have therefore the same entanglement properties. 

As an example, by considering a graph of four qubits linked in a row, the corresponding 
cluster is given by 

I*?"") = ^(1 + 00+) + I + 01-) + I - 10+) - I - 11-)) . (1) 

One way of generating cluster states is using photons. A useful tool to realize multiqubit 
states is represented by the so-called hyperentanglement (HE), i.e. the entanglement of two 
(or more) particles in several degrees of freedom (DOFs) [2T]. Precisely, by using the source 
shown in figure [l} we can generate two photons hyperentangled in polarization and path: 

IS4) = -^mjH)^ + \V)JV)^) ® ^(K) Jr)^ - |r) J£)^) . (2) 

In the previous equation \H) {\V)) represent the horizontal (vertical) polarization state and 
|r) and \i) are the two modes (right and left) in which each photon (^4 and B) can be emitted. 

The two photons (at degenerate wavelength A = 728 nm) are emitted by the spontaneous 
parametric down conversion (SPDC) process in a nonlinear type-I ;5-barium-borate (BBO) 
crystal. The BBO crystal is shined by a vertically polarized continuous wave (cw) Ar+ laser 
(A^ = 364 nm). Polarization entanglement is generated by the double passage (back and 
forth, after the refiection on a spherical mirror) of the UV beam. The backward emission 
generates the so called V-cone: the SPDC horizontally polarized photons passing twice through 
the quarter waveplate (QWP) are transformed into vertical polarized photons. The forward 
emission generates the H-cone. Temporal and spatial superposition are respectively guaranteed 
by the long coherence time of the UV beam and by aligning the crystal at a distance from the 
spherical mirror which is equal to its radius of curvature. In this way, the indistinguishability 
of the two perpendicularly polarized SPDC cones creates polarization entanglement: when two 
photons are detected it is impossible to know in which pump passage through the crystal they 
have been generated. It is worth noting that the probability of double pair emission (i.e. four 
photons) is negligible due to the low power of the cw pump beam (< 100 mW). By translating 
the spherical mirror it is possible to change the relative phase between the states \HH)ab and 
a lens L located at a focal distance from the crystal transforms the conical emission 
into a cylindrical one. 

Path entanglement can be generated by exploiting the properties of Type-I phase 
matching. The two polarization entangled photons are emitted over two opposite directions of 
the SPDC cone. By selecting with a four-holed mask two pairs of correlated modes, thanks to 
the spatial coherence property of the source, the two photons are also entangled in path. We 
labeled the two pairs of correlated modes, as — and ^a — '^b- We set the relative phase 
between the two pair emissions to the value 99 = tt by tilting think glass on the photon paths. 
The state expressed in ^ encodes 4 qubits into 2 photons [22l |23] . 
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2.1. 4-Qubit cluster 

The hyperentangled state may be transformed into a cluster state {C^) by using a waveplate 
with vertical optical axis and placed on the |r) mode of the A photon. The waveplate acts as 
a TT phase shift on the state \V)B\r)B- When applied to IS4) it generates the following cluster 
state: 

1^4) = lim^Hr)^ - \Hr)jm)^ + \Vi)jVr)^ + \Vr)jVi)^) . (3) 

By using the correspondence \H)a,b ^ |0)3,4, \V)a,b ^ |1)3,4, \£)a,b ^ |0)2,i, \r)A,B ^ |l)2,i, 
the generated state {C^} is equivalent up to single qubit unitaries to |$4^): {C^) = U\^^}^) = 
XiHi (g) Z2 I3 1 $4^), where H represents the Hadamard gate H = -^{X + Z). The 

latter relation between \C4) and |$4^^), implies that \C4) is the only common eigenstate of the 
generators % = UgaU~^ obtained from ga by changing Xi ^ Zi^ Zi ^ —Xi^ X2 —X2 and 
X4 ^ Z4. 

2.2. 6-qubit cluster 

It is possible to add more qubits to the state by selecting more optical paths. Precisely, by 
selecting four pairs of modes it is possible to generate a two-photon six-qubit hyperentangled 
state [24j. We labeled the four modes on which each photon can be emitted as \Er)^ 
\Ir) and |/^), where E (/) stands for external (internal) mode (see figure [2^)). The 6-qubit 
hyperentangled state can be written as 




Figure 2. Generation and measurement of the 6-qubit linear cluster, (a) By selecting 4 
pairs of correlated SPDC modes, a 6-qubit polarization-path hyperentangled state can be 
generated. Two half waveplates (A/2) are used to transform it into a 6-qubit linear cluster 
state corresponding to the graph shown in the inset, (b) Two cascade interferometers are 
used for path measurement. The first BSi performs the measurement in the qubit 
for both photons, while the two beam splitters BS2A and BS2B perform the measurement 
in the {|/), \E)} qubit for Alice and Bob photon, respectively. Each detection stage Di is 
composed of a polarization analyzer (waveplates and polarizing beam splitter) followed by a 
single photon detector. Two translation stages change the optical delay Axi^2 to obtain the 
correct temporal superposition of the different modes. 
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iSe) = —mjH)^ - \V)JV)^) ^(|r) + \i)Jr)^) 



®^{\E)^\E)^ + \I)JI)^). (4) 

As shown in figure [2^), two half-waveplates are used to transform the previous state into the 
6-qubit hnear cluster state [71 [12] : 

iLCe) = ^[\EE)\cl^^Uri) + \EE)\rUir) + |//)|^+).|r^) - (5) 

where \^^)^ = 1/^2(1 W) ± |ym and 10^)^ = 1/V2{\HH)±\VV). \LCg) corresponds to the 
graph shown in the inset of figure 2pL) up to single qubit unitaries. Precisely, iLCe) is the only 
common eigenstate (with +1 eigenvalues) of the generators gi obtained from gi by changing 
X2 ^ 2^2, X3 ^ — Z3, Z3 ^ X3, X4 ^ Z4 and X5 —X^. In order to measure Pauli path 
operators, two cascade interferometers are implemented (see figure [2): 



3. Results 



3.1. Quantitative Entanglement Verification 

The detection and quantification of entanglement has become a standard part of quantum 
information experiments. Methods for entanglement detection range from Bell inequalities 
over entanglement witnesses [25j to semidefinite programs [261 I27j . In order to quantify 
entanglement, it is necessary to evaluate an entanglement measure for the state under scrutiny 
[13j. Entanglement measures have the advantage that they do not only detect entanglement, 
but they may also provide an operational meaning to the amount of entanglement in a given 
quantum state. Until today, many entanglement measures have been invented, and the choice 
of the appropriate measure depends on the specific task [T3] . 

Here, we choose the global robustness of entanglement [T7l[l8l[l9] and the relative entropy 
of entanglement [20j . Both measures are suitable to quantify graph state entanglement for the 
following reasons: cluster states were introduced as multipartite entangled states that exhibit 
a particular persistence against noise. While GHZ states become more vulnerable under noise 
with increasing system size, this is not the case for cluster states [T]. Hence, the robustness 
of entanglement is a measure that quantifies this property. The relative entropy provides an 
operational meaning for cluster states in the sense that it 'counts' the number of entangling 
gates. As shown in Ref. [35l [28j the relative entropy of entanglement for cluster states is 
proportional to the number of applied cont rolled-phase gates. 

The relative entropy of entanglement is defined as j20] 

Er{p) = min S{p\a), (6) 
where Sep denotes the set of fully separable states, and S{p\a) = tr[p(log2 p — log2 cr)]. 
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Figure 3. Raw and optimized values of p^i^y for the 4-qubit cluster state. Each value 
corresponds to a given p^i^y according to the correspondence {k} ^ Sk^ + 4/^3 + 2/c2 + ki. 



The global robustness of entanglement measures how much noise must be mixed in to a 
given quantum state such that the mixture becomes separable [TTl [HI [I9] : 



7^ / \ • r p + sa ^ . 
RgKP) = mm {s : G bep\, 



(7) 



where V is the entire Hilbert space. 

From a mathematical point of view it is more convenient to relax the global robustness by 
replacing the set of fully separable states by teh set of PPT states, thus obtaining the following 
semidefinite program: 

= mintr{a} (8) 
subject to a > 0, (9) 

[p + (yf >o. (10) 

Here F denotes partial transpostion with respect to a partition of choice. In principal, one 
could check all possible partitions. In this way we have relaxed the global robustness to a 
PPT version that can be formulated as a semidefinite program. Hence, numerical tools such 
as convex optimization solvers are instantly available to evaluate this measure [29j. 



3.2. Four-quhit cluster state 

To verify the creation of the four-qubit cluster state all elements of the stabilizer group 
were measured [30j. As the measurements were local measurements, also statistics of single 
Pauli operators are available. These local measurements do not contribute to the fidelity, 
but they allow us to improve bounds on entanglement measures as those are restricted 
minimizations that can only improve when more constraints are added [32j. Using the 
measured data we calculated the raw fidelity F\c^) = jQj2{ka}=o('^^ik'2k3k4) = 0.880 ± 0.013 
[30] . where Skik2k3k4 = Ylt=ii9o')^'' stabilizers built as all possible products of 

generators. The raw purity is found to be P{p) = tr{p^) = 0.779 ± 0.005. From the raw 
data it is possible to obtain the fidelity with all possible graph state bases |{^})({^}| since 
|{^})({^}| = Te "^{k} '^{k}{~^y^ with i • k = "^21=1 ^aka- Some of the raw fidelities are negative 
because of experimental inaccuracies and statistical fluctuations of coincidence counts (the 
same problem arises in quantum state tomography \3T\). To solve the problem we applied a 
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Table 1. Four-qubit cluster state: measurement results of generators 

Generator Measurement Outcome 

gi = -Zi (g) Z2 (g) I3 (8) I4 0.994 ± 0.001 

g2 = -Xi (g) X2 (g) Z3 (g) I4 0.849 ± 0.003 

^3 = li g) Z2 g) X3 g) X4 0.937 ± 0.003 

^4 = li g) I2 g) ^3 g) Z4 0.911 ± 0.002 



maximum likelihood estimation. We determine the physical density matrix diagonal in the 
graph state basis and written as 

pphys = X]ml{^})({^}l ' m ^ (11) 

{k} 

that is most compatible with the experimental data. The value of the optimized poooo 
corresponds to the fidelity with the cluster state, with value equal to 0.880 which is completely 
compatible with the raw fidelity. The other values of the optimized p^k} are shown in figure 
jS, The optimized purity is given by P^pt = J2{k}P{k} ^ 0.778 again compatible with the raw 
value. 

While for few-qubit systems the determination of the whole stabilizer is feasible, this is 
not the case for large graph states. Therefore, it is natural to ask which bounds on the fidelity, 
purity and entanglement can be obtained from incomplete information on the density matrix 
[321 [33I [34] . e.g., from measuring generators of a graph state only; and how do such bounds 
scale with system size under realistic conditions? 

Let us consider the estimation of the fidelity from information on the generators only. 
This is formulated as a worst-case estimation [14j 

Fmin = min{F(p) : tr{pgi) =ai,p> 0}, (12) 
p 

where g^ are the generators of the graph for i = 1,..,4 with corresponding measurement 
outcomes a^, and go = 1. Remarkably, this problem can be solved optimally, leading to a 
solution of the analytic form 



F_ = max{0, SkhLji±^} (13) 

for n qubits and holds for all stabilizer operators with spectrum {+1,-1}. One may quickly 
check that the optimal lower bound on the fidelity consistent with the measurements of the 
generators (see Tab. ([T])) is given by F^in = 0.846 ± 0.009. The relative loss of information 
on the fidelity is therefore only around 5%, even though only four out of sixteen elements of 
the stabilizer were determined. It is also possible to optimally estimate the purity using only 
generator measurements. Following the techniques of Ref. [15] we obtain a minimal purity 
consistent with such measurements of Pmin = 0.715±0.014. 

Quantifying the experimentally created entanglement is achieved by evaluating the global 
robustness of entanglement and the relative entropy. The density matrix reconstructed from 
the stabilizer measurements and local observables obtained as a side product of the stabilizer 
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measurements serves as the input for the semidefinite program ([8|), where we evaluate the 
global robustness with the constraint of positivity of the partial transpose with respect to all 
partitions. We find that the PPT-Robustness is given by = 2.519 ± 0.012. Note that 

this value represents a lower bound to the global robustness in its standard version ([7|. It is 
straightforward to compute the logartihmic global robustness: LRq^^ = 1.817 ±0.005, which 
is not far from its desired value of 2. 

As in the case of the fidelity one might ask which bound on the entanglement can be 
obtained from generator measurements alone. Here the estimation is analogously formulated 
as a minimization of the measure over states consistent with the measurement data 

= min{i?G(p) : tr{pg,) = a,, p > 0}. (14) 
p 

A lower bound to this problem was derived in Ref. [T6], namely Rcmin = 
max{0, 2l^l(^^^^-^^^^) — 1}. Here B denotes the smaller set of qubits resulting from a 
coloring of the system into two colors, say Amber A and Blue B with \A\ > \B\ (see Ref. [33] 
for more details). With this formula, we attain the following analytic bound on the global 
robustness based on the outcomes of the generators only [16]: 

= 2.384 ± 0.036. (15) 

In turn, one can then easily compute an analytic bound on the logarithmic global robustness: 
LRg = 1.759. 

The relative entropy can also be bounded from below using techniques presented in Ref. 
[IB]. The problem reads: 

ERrmn = mm{ER{p) : tr{pg,) =a,,p> 0}. (16) 
p 

A lower bound to this minimization is given by 

E^_=max{0,|S|-^//(p,)}, (17) 

i 

where pi = and H{x) = — xlog(x) — (1 — x) log(l — x) is the classical entropy function. 
Then, by merit of Eq. ( [T7| ) we achieve the following bound on the relative entropy of 
ent an dement: Eu = 1.120 ± 0.021. Using; all stabilizer measurements we find a lower 
bound to the relative entropy of 1.449±0.013. Hence, the relative difference of the entanglement 
bounds of the relative entropy is considerably larger than the relative difference of the estimate 
of the robustness to its real value. 



3.3. Six-qubit cluster state 

The six-qubit cluster state is verified utilizing the same techniques as in the four-qubit case. 
All 64 stabilizer operators were measured and mapped to a density matrix via maximum 
likelihood, giving a fidelity of F = 0.645 ± 0.006. Estimating the fidelity from the generators 
alone gives Fmin = 0.545 ± 0.027. 

Next, we obtain for the purity P{p) = tr{p^) = 0.424 ± 0.010, and the worst-case purity 
estimate from the generators Pmin = 0.297 ± 0.015. We compute the PPT-Robustness for 
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Generator Measurement Outcome 



91 


= Xi (g) X2 (g) I3 


C 


^15 (^le 


0.593 ±0.008 


92 


= Zi (g) Z2 (g) I3 g) 


I4 g) 


Z5g)l6 


0.879 ± 0.005 


93 


= -li g)l2 g)Z3 


g) I4 ' 


g) I5 (8) Ze 


0.998 ±0.001 


94 


= Zi g) I2 g) I3 g) 


Z4 g) 


I5 g)l6 


0.997 ±0.001 


95 


= -li g)X2 g)l3 


g) I4 


g) X5 g) Ze 


0.791 ±0.006 


96 


= ll g)l2 g)X3 (g 


) I4 g) 


Z5 Xq 


0.831 ±0.006 



the reconstructed state and find Rq^^ = 4.507 ± 0.047 resulting in a logarithmic PPT- 
Robustness of LRq^^ = 2.461. If one only measures the generators of the stabilizer, one 
obtains Rcmin = 3.360 ±0.216 and LRq^^^ = log(l + = 2.124 respectively. This means 

that despite the lower fidelity and lower fidelity estimate one obtains a higher bound on the 
entanglement, even though only 6 out of the 64 elements of the stabilizer are used to obtain 
the bound. 

To obtain a bound on a second measure, the relative entropy of entanglement, we use Eq. 
([TtI ) to obtain Er^.^ = 1.013 ± 0.046. Using all stabilizer measurements, one obtains a lower 
bound to the relative entropy of 1.492 ± 0.027. 

4. Conclusion 

We have presented the creation of four- and six-qubit cluster states using photons. The 
cluster state entanglement was encoded in path and polarization DOF, thus rendering the 
state hyper entangled. A summary of the relevant characteristics of the created states is given 
in Tab. [3| The created state could serve as the basis for one-way quantum computation and 
represents an important step in realizing optical quantum computing. 

The experimentally created entanglement was quantified in terms of entanglement 
measures, namely the global robustness of entanglement and the relative entropy of 
entanglement. Our results also give insight into the question how analytic bounds from 
incomplete tomographic information scale with the system size under realistic noisy conditions. 
Our results demonstrate that despite the decreasing fidelity and purity of the state, one can 
still infer higher amounts of entanglement with a number of observables which is linear in the 
number of constituents. 



Table 3. Summary of experimental results 





Fidelity 


Purity 


Minimal Global Robustness 


Minimal Relative Entropy 


4 qubits 


0.880 ± 0.006 


0.778 ± 0.005 


2.519 ±0.012 


1.449 ±0.013 


6 qubits 


0.645 ± 0.006 


0.424 ±0.010 


4.507 ±0.047 


1.492 ±0.027 
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